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ABSTRACT 



We present initial result of a large spectroscopic survey aimed at measuring the timescale of mass accretion in young, pre-main- 
sequence stars in the spectral type range KO - M5. Using multi-object spectroscopy with VIMOS at the VLT we identified the fraction 
of accreting stars in a number of young stellar clusters and associations of ages between 1-50 Myr. The fraction of accreting stars 
decreases from ~ 60 % at 1.5 - 2 Myr to ~ 2 % at 10 Myr. No accreting stars are found after 10 Myr at a sensitivity limit of 10~" 
M^yr' 1 . We compared the fraction of stars showing ongoing accretion (face) to the fraction of stars with near-to-mid infrared excess 
(fiRAc)- In most cases we find f acc < f IRA c > i.e., mass accretion appears to cease (or drop below detectable level) earlier than the dust is 
dissipated in the inner disk. At 5 Myr, 95 % of the stellar population has stopped accreting material at a rate of > 1 0~ M Q yr~ 1 , while 
~ 20 % of the stars show near-infrared excess emission. Assuming an exponential decay, we measure a mass accretion timescale (t„ cc ) 
of 2.3 Myr, compared to a near-to-mid infrared excess timescale (t IRA c) of 2.9 Myr. Planet formation, and/or migration, in the inner 
disk might be a viable mechanism to halt further accretion onto the central star on such a short timescale. 

Key words. Accretion disks - Stars : pre-main-sequence - Planetary systems: protoplanetary disks 
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1. Introduction 

Circumstellar disks around young, pre- main-sequence stars are 
the natural loci of planet formation (e.g. lHennindl2008l) . Such a 
protoplanetary disk is formed during the collapse of the molec- 
ular cloud to conserve the initial angular momentum. It is made 
of interstellar gas and dust. The conventional pla net formation 



model is the so-called core - accretion model (e.g. Pollack et al 



1996: ]Mordasini et al.ll2008l) . In the simulations of[Pollack et al 
(1996), the formation of solar-system like planets takes up to 16 
Myr. Giant planet form ation is much faster (~ 3 Myr) if planet 
migration is included ( Alib ert et alj|2004l 120051) . Gravitational 
instability in the disk was also proposed as a viable scenario 
to form plane ts, especially in the outer disk dBosslll997l Il998t 
Durisen 2007). If gravitational instabilities occur, the disk may 
fragment and form dense self-gravitating clumps which are the 
precursor of gas giant planets. In this model planet formation is 
very fast (< 1 Myr). 

Planet formation models face the fast evolution of protoplan- 
etary disks. Star forming regions show a steady decline with 
time in the fraction of stars having infrared excess emission (e.g. 
Haisc h et al. 1 1200 1L 12005b iouwmanitiD 120061: iHillenbrand 



2008). This is thought to be caused by gradual clearing of dust 



in the inner disk. The growing consensus is that the warm small 
dust grains disappear within the first 5-^9 Myr (Tdust)- The quan- 
tity Tdust is sometimes adopted as the disk lifetime. Although 
dust is essential to form a planet, in terms of disk mass, and 
hence dynamics and evolution, it accounts only for a small per- 
centage. The bulk of the disk mass is thought to be gaseous. At 



present, it is unclear whether the gas dissipation timescale (T gas ) 
is similar to Tdust- 

In this paper we report on a study aimed at determining the 
timescale for mass accretion (T nee ) in protoplanetary disks. The 
timescale r acc , i.e. the time at which the disk accretion phase 
ceases, provides a strong constraint on T gas . Gas may still be 
present after T acc , however, the amount of remaining gas, and 
hence of disk mass, may be too low to be able to form giant 
planet(s). Our observational strategy, observations and data re- 
duction are presented in section 2. Analysis and results are pre- 
sented in section 3 and 4 respectively. In section 5 we compare 
our result with literature data. Conclusions are drawn in section 
6. 



2. Method 

In order to measure r acc in a secure way, we performed an opti- 
cal spectroscopic survey of a large number of stars towards seven 
young stellar clusters (Table [TJ. Multi-object s pectroscopy was 
performed with VIMOS3 dLeFevre et alj|2003l) at the VLT. No 
a priori target selection was made, but instead we optimized the 
position of the slits in the multi-object masks in order to get the 
spectrum of as many stars as possible in the proximity of the 
cluster center. We did not aim at taking the spectrum of each 
(known) cluster member, but rather at building an unbiased in- 
ventory of the stellar population in a sub-region of the cluster 
down to a limiting magnitude of V = 21 mag. For each cluster 
we estimated the fraction of stars with ongoing mass accretion 
by analyzing the Ha line profile. Finally, the fraction of accret- 
ing stars is plotted against the age of the clusters. In this paper 



* Based on observations collected at the European Southern 
Observatory, Paranal, Chile (Proposal ID: 078.C-0282; 081.C-0208) 



http://www.eso.org/sci/facilities/paranal/instruments/vimos/ 
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Table 1. Regions observed with VIMOS. The surveyed area for each cluster is 16'xl8' centered at the coordinates listed below. The 
limiting magnitude is V = 21 mag. MOS observations were carried out with the "HR-Orange" grism (maximum spectral coverage: 
5200 - 7600 A; A/ AA = 2 1 50) and a slit width of 1 " . 



Cluster 


Observing 


period 


RA(J2000) 


DEC(J2000) 


Distance 


A v 


Sources 








[hh:mm:ss] 


[dd:mm:ss] 


[pc] 


[mag] 


[N] 


a Orionis 


Oct. 2006 


- Mai'. 2007 


05:38:00 


-02:37:00 


360 


0.15 


216 


NGC 6231 


Apr. 2008 


- Sep. 2008 


16:54:10 


-41:49:30 


1240 


1.36 


573 


NGC 6531 


Apr. 2008 


- Sep. 2008 


18:04:13 


-22:29:24 


1210 


0.87 


627 


ASCC 58 


Oct. 2006 


- Mar. 2007 


10:15:07 


-54:58:12 


600 


0.28 


370 


NGC 2353 


Oct. 2006 


- Mar. 2007 


07:14:30 


-10:16:00 


1120 


0.22 


316 


Collinder 65 


Oct. 2006 


- Mar. 2007 


05:26:05 


+ 15:41:59 


310 


0.40 


198 


NGC 6664 


Apr. 2008 


- Sep. 2008 


18:36:42 


-08:13:00 


1440 


2.20 


508 



we present the analysis of the mass accretion evolution. Results 
on individual clusters will be presented in a forthcoming paper. 
The following criteria were adopted for the sample selection: 

Age : the clusters have ages covering the interval 2-30 Myr; 
Distance and extinction : nearby and optically visible systems 

were selected to ensure detection of low mass stars down to 

spectral type M3 (at least); 
Stellar density profile : to avoid any bias introduced by, e.g., 

mass segregation, we only selected clusters known to have 

a spherical density profile; 
Angular radius : (also linked to previous point) only clusters 

with angular size comparable to the field of view (FOV) of 

VIMOS were considered in order to minimize the number of 

contaminants. 

For each cluster an area of 16'xl8' was investigated acquir- 
ing spectra of ~ 200 - 600 representative objects in the central 
region of the cluster. Multiple dithered exposures were taken. 
Details about observations are given in TableQ] 

Data reduction was performed with IRAF and specific IDL 
routines. The spectra were reduced using a long-slit approach, 
i.e., extracting each spectrum separately. Due to distortion of 
the instrument night time wavelength calibration exposures were 
taken. No night flat field images were taken. For this reason we 
preferred to not correct the MOS spectra for flat field to not in- 
troduce further source of noise. The following scheme was ap- 
plied: bias subtraction, wavelength calibration, background re- 
moval and s pectrum extr action using the Home optimal extrac- 
tion method dHornelll986l) . Multiple dithered spectra were com- 
bine afterwards using a specific IDL routine: before combining 
the spectra, a scaling factor was applied to take into account slit 
losses in different exposures. 

Fig.Q]shows some examples of VIMOS spectra of our pro- 
gram stars. 

3. Analysis 

3. 1 . Spectral type 

The spectral type of the objects classified as cluster members 
was derived. For this purpose we used the TiO index defined as: 

CY7020 - 7050A) 

770(7 140A) = — '— (1) 

770(7125 -7 155 A) 

where C(7020 - 7050 A) is the spectral continuum com- 
puted between 7020 - 7050 A and TiO(7125 - 7155 A) is the 
intensity of the TiO molecular band absorption between 7125 - 



7155 A (lAllen & StrorrJl995l:lBricefio et al.ll998tlOIiveira et alj 
120031) . The TiO molecule is very sensitive to the stellar gravity 
and correlates very well with the spectral type in the range K5 
- M6. We computed the spectral type of th e VIMOS spectra, 
following the scheme of iJeffries et al.l (120071) . For earlier spec- 
tral types (K0 - K5) we directly compared the VIMOS spectra 
with a set of standard stars spectra covering the same spectral 
range and having similar spectral resolution. The uncertainty on 
spectral type derived with these techniques is one sub-class. 

In the following analysis we consider only stars having spec- 
tral type between K0 and M5. 

3.2. Membership 

The cluster membership was established on the basis of the pres- 
ence of two spectral features: emission of the Ha line at 6563 
A and absorption of the Li 6708 A line. Ha emission in young 
stars is either produced by mass accretion or by chromospheric 
activity. The presence of Li is also widely used as youth indica- 
tor (e.g. iBasri et al.l[l99ll) : Li is depleted very fast once the ther- 
monuclear reactions begin in the core of the star. The abundance 
of Li in the stellar atmosphere i s inversely proportion al to the age 
of the star ( IJeffries et al.ll2007h iMentuch et al.ll2008l) . This youth 
indicator is best suited for K- and M-type stars younger than 
~ 30 Myr (IJeffries et alj|2007l) . Some stars have strong Li ab- 
sorption in their spectrum but no Ha emission. These might be 
cluster members with no, or reduced, chromospheric activity. We 
measured the EW of Li 6708 A of these sources and compared 
this with the typical Li 6708 A EW measured for the certain 
members of the cluster (0.2 A < EW [Li] < 0.6 A). We note 
that the lower limit of EW Li measured here (0.2 A) is an inter- 
mediate value assumed as Li indicator applied by other authors 
(e.g. IJeffries et al.ll2007l EW > 0.3 A: IWinston et~ai1l2009L EW 
> 0.1 A). In most of the cases the spectra show either a tiny (EW 
< 0.1 A) Ha emission or small absorption. In the latter case, 
the absorption is much smaller than the typical photospheric ab- 
sorption for the same spectral type. The "reduced" Ha absorp- 
tion might result from line veiling caused by the chromospheric 
activity. These stars are classified as cluster members. 

3.3. Signature of mass accretion 

To distinguish between chr omospheric an d accre tion ori- 
gin of the Ha emission, IWhite & BasrH (|2003) suggest 
a different threshold in EW for different spectral types. 
iBarrado v Navascues & Martini ([2003 ) have defined an empiri- 
cal criterion based on the saturation limit of the chromospheric 
activity. Another tool of investigation is the width of Ha at 
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10% ( Hai 0% ) of the line peak (|^Qiite_&Bas ri 2003; Natt a et alj 
I2004t Uavawardhana et all l2006t iFlahertv & Muzerolld [2Q08h . 
As shown in Fig. 7 of IWhite & Basril d2003l) . accreting objects 
have systematically larger value of Haio% compared to non- 
accreting objects. Stars with Haio% > 270 kms 1 are accreting 
material from a circumstellar disk. This criterion is valid down 
to a m ass accretion rate of 10 -10 4- 10 -12 M yr _1 dNatta et al.1 
120041) . The actual limit on detectable accretion rate depends on 
the stellar mass and age. As a conservative value we adopt 10~ u 
M Q yr~ l . In this work we used both EW [Ho-] and Haio% cri- 
terion to distinguish between accretion and chromospheric ac- 
tivity. This method is solid to find signature of accretion when- 
ever EW \Ha] and Hotiq^ are above, or close to the thresholds 
dJavawardhana et all 120061) . However, stars with values lower 
than the cutoff values may still be accreting at low rate (< 10~ u 
M )T -1 ) and cannot be distinguished from pure chromospheric 
activity. 

The spectral resolution of VIMOS is enough to measure the 
widening of the Ha line due to accetion. At 6563 A the minimum 
velocity resolvable with VIMOS is 140 kms _I . We measured 
the EW and the 10% width of the Ha emission. The spectral 
continuum was computed using two spectral windows at both 
sides of the Ha line. The continuum level in correspondence of 
the emission line was interpolated using a first order polynomial. 
The spectra were then continuum-normalized and subtracted. 
Haio% was measured as the width of the Ha line at the height 
corresponding to 10% of the line peak. In the case of narrow 
emission line and low signal-to-noise spectrum, the 10% level 
might be confused with the neighbor continuum fluctuations. To 
avoid this, the result was visually inspected. The major source of 
uncertainty for the measurements of EW [Ha] and Haio% is the 
determination of the underlying spectral continuum. The errors 
on Ha EW and Haio% were computed by taking the continuum 
error into account. For most of the spectra the signal-to-noise ra- 
tio close the Ha line is > 10 - 20. This translate in a precision 
of ~ 0.1 A and 5 kms -1 for equivalent width and Haio%, re- 
spectively. A further contribution of uncertainty is given by the 
instrument's spectral precision. We determined this by comput- 
ing the full widht half maximum and standard deviation for a 
number of arc lines in proximity of the Ha line. The standard 
deviation is of the order of 5 - 7kms~ I . We therefore assume 
0.2 A and 10 km s _1 as a (conservative) lower limit to the uncer- 
tainty of EW [Ha] and Ha 10 %, respectively. 



4. Result 

None of the stars in regions older than 10 Myr in our survey 
(ASCC 58, Collinder 65, NGC 2353 and NGC 6664) shows ev- 
idence of ongoing accretion (Table |2}. Results for the youngest 
(< 10 Myr) regions are described below. In Fig. [2] we plot EW 
[Ha] versus Haio% for three different spectral ranges. Vertical 
and horizontal dashed lines represent the thresholds for accret- 
ing stars. The different symbols refer to different clusters : dots 
- cr Ori, triangles - NGC 6531, stars - NGC 6231. 

As shown in Fig. [2] there are a number of sources with 
EW [Ha] below the accretion threshold but with Haio% > 270 
kms -1 . If the line broadening is caused by mass accretion at a 
rate > 10~ u M Q yr~ l , we would expect to measure a large EW as 
well. Beside accretion, stellar binarity and fast rotation might 
also be responsible for the line broadening. In order to deci- 
pher the nature of these sources we also investigated the pres- 
ence of other emission lines which are associated with mass 
accretion (e.g. Hei 5876A, 6678A; e.g. iBeristain et all l200lt 



iHerczeg & Hillenbrandl20 08). In the presence of these accretion 
diagnostics, we can safely state that the star is accreting. 



o-Ori 

10 out of the 216 sources with spectra taken with VIMOS are 
classified as members of the crOri cluster. All of them show 
emission in the Ha line as well as presence of Li. 2 out of 10 
objects are clearly accreting stars according to the method men- 
tioned above. A further object (<x Ori # 1: RA = 84.584154, 
DEC = -2.633792, Sp. Type M5) has a broad Ha 10 % (350 
km s" 1 ) but an Ha EW value (15 A) below the accreting thresh- 
old. The presence of He remission lines ( 5876 A EW = -0.75, 
6678 A EW = -0.2) suggests the presence of hot gas in the 
vicinity of the star a nd hence signature of ongoing accretion 
(IBeris tain et al. 2001). The fraction of accreting objects in <x 
Ori derived here is 3/1 or 30 % + 17 %fl Th is is in agreement 
with previous result of IZapatero Osorio et al.l (|2002|, 30-40 %), 
iBarrado v Navascues & Martini d2003l 27 + 7 %). lOliveira et al.l 
(2006, 30 %). The good agreement with previous results vali- 
dates our observational strategy despite the low number statis- 
tics. 



NGC6231 

Out of 573 sources in the region of NGC 6231 we identified 78 
objects as cluster members. The majority of these sources have 
spectral type in the range K0 - M5 (there are 3 sources of earlier 
spectral type) and show both Ha emission and strong Li absorp- 
tion. 1 1 out of 78 objects (all late K type) were identified based 
only on the presence of Li. Each of these spectra was compared 
with a standard star of equal spectral type. All these spectra show 
either a tiny (EW ~ 0. 1 A) Ha emission or small absoption. In 
the latter case, the absoption is much smaller than the typical 
photospheric absorption for the same spectral type. Further 3 
objects have no Ha and the presence of Li 6707.8 absorption 
is not clear. 7 out of the 75 K0 - M5 sources in NGC 623 1 are 
consistent with ongoing mass accretion with both diagnostics. 
Further 2 sources have EW [Ha] close to the accretion threshold 
but large (> 270 km s" 1 ) Hai % (panel K0 - K5 in Fig.|2|. Many 
stars with spectral type in the range K6 - M2.5 have large Haio% 
but small EW [Ha]. We inspected all of them and found 2 objects 
with spectral type K7 showing Hei 5876A in emission (EW = - 
0.5A, -0.6A respectively). The evidence of large Haio% together 
with the He i emission is most likely due to ongoing mass ac- 
cretion and these two stars are classified as accreting stars. We 
estimate a fraction of accreting stars in NGC 6231 of 11/75 or 
15 (+ 5 %). We warn the reader that this might be a lower limit 
to the actual fraction of accreting stars; further investigation is 
needed to disentagle the nature (accretion vs binarity/rapid rota- 
tion) of the systems with large Haio% (> 300 kms -1 ) but low 
EW [Ha]. 



NGC6531 

We identified 26 cluster members in NGC 6531 based on the 
presence of Ha emission and presence of Li. Further 13 sources 
show presence of Li 6708 A, but have Ha in absorption. As 
in the case of NGC 6231, these might be cluster members with 
no, or reduced, chromospheric activity. We measured the EW of 
Li 6708 A of these 13 sources and compared with the typical 



computed as JN, t accIN,, 



D. Fedele et al.: Accretion Timescale in PMS stars 




a 

3 
O 

o 



SpT : M4 

EW[Ha] : -15A 

Ha. n ^ : 350 km s" 1 
lu%> 



SpT : K4 
EW[Ha] : -33A 



'I '"'" ' "' 1 1 1 1 



SpT : M4 

EW[Ha] : -4.8A 

Ha iriw : 225 km s"' 
10% 



SpT:K5 

EW[Ha] : -7 A 

Ha. n/W : 700 km s" 1 
1U% 



SpT : M2 
EW[Ha] : -4.9A 




SpT : K7 

EW[Ha] : -2.6A 

~Ra. nrrf : 300 km s" 1 
1U% 



SpT : Ml 
EW[Ha] : -2.3A 



SpT : M4 

EW[Ha] : -4.8 A 

Hfl. nw : 230 km s" 
10% 




SpT : Ml 

EW[Ha] : -25A 

Ha 1 , w : 620 km s" 
10% 
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X [Angstroms] 

Fig. 1. Some examples of VIMOS spectra 



EW of the 26 stars in NGC 6531 showing also Ha emission (0.2 
A < EW Li < 0.6 A). 10 of the 13 sources have similar EW 
and are likely cluster members. The remaining 3 sources have 
EW between 0.1- 0.2 A. Moreover, these are characterized by 
strong absorption lines at 5778, 5796, 6284, 6614 A which are 
produced by diffuse interstellar absorption bands (DIBs). These 
are likely located in the background of the cluster and are clas- 
sified as non-members. Out of the 36 identified members, 3 ob- 
jects show broad and strong Ho- emission consistent with mass 
accretion. One further object is classified as uncertain accreting 
star. The fraction of accreting stars is then 8+5%. Given the 
large Haio% (> 300 kms -1 ) of 2 objects in the spectral range 
M3 - M5, the fraction of accreting stars computed here might be 
a lower limit. Follow up observations are needed for these two 
sources. 



5. Discussion 

The fraction of accreting stars (face) f° r eacn cluster is listed in 
Tableland plotted in Fig.[3]together with literature data. 

The aim of this paper is to trace the evolution of mass ac- 
cretion with time. A critical point is the age of the clusters. We 
adopt ed the most recen t age measurements (Table For NGC 
623 1 ISana et al.l d2007l) find that the bulk of the CTTS has an age 
between 2-4 Myr, although with a large age spread. The most 
recen t study of cr Ori suggests an age o f 3 Myr (e.g. ICaba llero 
120081) . For NGC 6531 iPark etail d2001l) me asure an age of 7.5 
+ 2 Myr in agreement with (Forbes 1996) who measure 8 + 
3 Myr. For the remaining clusters we adopt the age provided 
by WEBDA0. For these clusters the age is estimated by fitting 



3 http://www.univie.ac.at/webda/, operated at the Institute for 
Astronomy of the University of Vienna 
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Fig. 2. EW [Ha] vs Haio% for three different spectral type ranges. Vertical and horizontal dashed lines represent the thresholds 
for accreting stars. The different symbols refer to different clusters : dots - Sigma Ori, triangles - NGC 6531, stars - NGC 6231. 
Objects classified as accreting stars are over plotted with open circles. 



the zero-age-main-sequence (ZAMS) to the brightest stars of the 
cluster. We take an accuracy of 30% as a conservative assump- 
tion. 

In Fig. [3] the fraction of accreting stars as function of clus- 
ter age is shown. The results obtained with VIMOS are shown 
as filled circles. Measurements of f acc exist in literature for 
some clusters. These are are shown as filled squares in Fig. 
[3] and are listed in Table [2] We considered only det ermina- 
tions of fa CC , using the method of IWhite & Basril d2003l) and/or 
Barrado y Navascues & Martin (2003). The fraction of accret- 



ing stars may vary w ith spectral type. In the case of Upper Sco 
Moh antv etaT] d2005l) e.g., found an increase from 7% (± 2%) 
for spectral type earlier than M4 to 20% (± 10%) for later spec- 
tral t ypes. For the younges t regions, p Oph, Taurus, Cha I and IC 
348, Moh antv et al.l d2005l) find no appreciable variation, within 
errors, in the fraction of accreting sources between the two mass 
ranges. For consistency with our survey we include here only the 
results for spectral types between K0 - M5. 



5.1. Evolution of f acc 

The fraction of accreting stars decreases quickly with time 
within the first 10 Myr. There is a clear trend from the 1.5 - 
2 Myr old regions (f acc = 60%) down to the 10 Myr old clus- 
ters (mean f acc = 2%). At the age of 5 Myr the mean f acc is 
~ 5% drastically lower than the 3^1 Myr old regions of cr Ori, 
NGC 6231 and Trumpler 37 (average 28%). No accreting stars 
are found beyond 10 Myr. This is in good agree ment with previ- 
ous measur ements of mass accre tion timesc ale by Mohantv et alj 
(e. g. 120051) : I Jav awai'dha na et al] (e. g. 120061) . There are some out- 
liers. The 77 Cha cluster shows an higher frequency of accre- 
tors c ompared to clusters of equal age (e.g. Sicil ia-Aguilar et all 
2009). This particular association shows a similar behavior in the 
disk frequency (from near-infrared excess). Age a nd distance are 
unlikely to be wrong for this well know system. iMoraux et alj 



(2007) suggest that most of the low-mass members of the as- 
sociation have been dispersed by dynamical evolution. The cur- 
rent list of members of r\ Cha might be biased towards infrared- 
excess/Ha emitting sources. Inversely the very young region 
of p Oph lies below the trend. The statistics in this region are 
limited by exti nction which may ob scure part of the stellar 



population (e.g . Mohantv et al.l [2005b . Moreover, the result of 
Mohantv et alTd2005l) might be contaminated by an older pop- 
ulation of stars; objects in the peri phery of p Oph core appear 
older than stars in the core itself ( Wilkin g et alj|2005h . In the 
following analysis we will not include p Oph and 77 Cha. 

NGC 623 1 (2 - 4 Myr) also appears to have a low f acc for its 
age. Interestingly, this region hosts many O and B stars. External 
photoevaporation (by means of the O, B stars) may accelerate 
disk dissipation. 



5.2. f acc vs fiRAc 

In this section we compare the fraction of accreting stars to the 
fraction of stars with near-to-mid infrared excess {fiRAc)- The 
latter is measured as the fraction of stars having infrared excess 
in the Spitzer/IRAC [3.6] - [8.0] bands over the total fraction of 
stars in the cluster. The definition of an excess source is based 
on the slope (a) of the infrared spectral ene rgy distribution d e- 
termined from the IRAC [3.6] - [8.0] bands (lLada et al.ll2006l) : 



oirac - 



d(log(AF A )) _ 
d(log(A)) ' 



3.6jim < A < 8.0pm 



(2) 



Star s with —1.8 < a < are classified as class 
II (e.g. iHernandez et al.l 120071) . In order to coherently com- 
pare fiRAc with face we refers only to stars of spectral 
type between K0 - M5. Whenever /irac was measured as 
function of stellar mass rather than spectral type we only 
consider stars of mass between 0.1 - 1.2 M Q . By apply- 
ing these selection criteria we find the following results for 
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Table 2. Adopted age, spectra l type range, f ac c and Jirac ( when available) in Fig. [3. 
S82). iParketalJ (f200Tl P OD. iHartmann et all d2005l Ha 05). iKharchenko etail (f2 005. K05), Moh antv et af] 
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Fig. 3. Accreting stars frequency as a function of age. New data 
(based on the VIMOS survey) are shown as (red) dots, litera- 
ture data as (green) squares. Colored version is available in the 
electronic form. 



- 50% dSicilia-Aguilaretail l2009h and NGC 7160 - 4% 
dSicilia-Aguilar et al.ll2006l) . 

Interestingly for some clusters we find a lower fractional 
value of stars with evidence of ongoing accretion compared to 
stars with near-to-mid infrared ex cess. For example am o ng the 
64 K0 - M5 stars identified by iFlahertv & Muzerolld (2008) 
in NGC 2068/71, 39 (61%) of these show clear signature of 
mass accretion while 45 (70%) stars have IRAC excess. We note 
that 2 of the non accreting stars in their sample are identified 
by st rong IRAC excess (ID 416, 843 in IFlahertv & Muzerollel 
2008), while one transitional object (ID 281) is a strong accre- 
tor. Similarly, for cr Ori we measure an f acc = 30% while fmAc — 
35%. This is consistent with the result of Damjan ov et al.l 12007) 



in Chameleon I. They find a small population of non-accreting 
objects bearing a dusty inner disk. 

In Fig. @] we plot f acc vs /irac for the regions where both 
quantities are known. Assuming an exponential decay we fit 
face and fmAc with the following function: 



fi = C ■ exp(-t I n) 



(3) 



where C is a constant and it is normalized assuming f acc = 
fiRAc - 100% for t — 0. The fit gives an accretion timescale of 
T acc = 2.3 Myr and a near-to-mid infrared excess timescale of 
t irac = 2.9 Myr. 



fmAc (Table 01: Taurus - 62% dHartmann et all l2005h . NGC 
2068/ 71 - 70% dFlahertv & Muzerollel l2008h. Cha I - 52% j 
64% dLuhmanet all 120081). IC 348 - 4 7% dLada et all 12006). 
cr Ori - 35% dHemandez et alj [2007b . Trum pler 37 - 47% 
dSicilia-Aguilar et alj|2006l), Upper Sco - 19% dCarpenter et alj 
2006), NGC 2362 - 19% dDahm & H illenbrandfl2007l). 77 Cha 



5.3. Inner disk dissipation and planet formation 

Mass accretion and dust dispersion in the inner disk appear to 
have different timescales. At an age of 5 Myr, 95% of the to- 
tal stellar population has stopped accreting material at a rate 
> 10~ n M yr~ l while ~ 20% of the disks still retain enough 
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Fig. 4. f acc (dots) versus /irac (squares) and exponential fit (dot- 
ted line for f acc , dashed line for /irac- 

dust to produce detectable infrared excess. In order to sustain 
mass accretion at a rate larger than 10~ n M & yr~ l a gas disk 
reservoir is required. This implies that at an age of r flcc = 2.3 
Myr most of the disk mass is drastically reduced while at an age 
of 5 Myr only 5% of the disks have enough gas mass to keep 
on accreting material onto the centr al star. Models of viscous 
evolution (e.g. Hartma nn et al.|[l998h predict a steady decrease 
of the mass accretion rate with time and this might be an expla- 
nation to the bimodal evolution of f acc and /irac- However, the 
timescale for viscous evolution to stop accretion (> 10 Myr) is 
much longer than our estimate of T acc (2.3 Myr). The efficiency 
of viscous evolution might be enhanced if photoevaporation i s 
also taken into account as first argued by Clar ke etall ((2001). 
Their so called "UV-switch" model is very efficient in stopping 
further accretion as the accretion rate drops below the photoe- 
vaporation rate. The advantage of this is the shorter accretion 
timescale. But this model also has a drawback, namely that disk 
dissipation following the end of the accretion phase should occur 
on a viscous timescale (10 5 yr). This is in contrast with the rel- 
atively high number of stars still retaining their dusty disk (e.g. 
20% at 5 Myr) after accretion ceases. 

Another possible explanation is that accretion onto the cen- 
tral star is stopped by planet formation, and migration, in the 
inner disk. If a giant planet forms/migrates to the inner region of 
the disk, it may stop further accretion onto the star. In this regard, 
recent s elf-consistent numeric al simulations of planet formation 
by, e.g jThommes et al" (2008) are consistent with this scenario. 
In these simulations planetesimal formation occurs over a large 
range of disk radii. As soon as the cores of giant planets form 
they migrate to small radii (< 1 AU). This happens very fast, 
within 1-2 Myr. The disk surface density drops to zero, and so 
does accretion onto the star, within the planet's orbit. Dust grains 
beyond the orbit of the planet instead, might still absorb and re- 
emit the stellar light and produce the measured IRAC excess. 

6. Conclusion 

In this paper we presented an analysis of the evolution of mass 
accretion in PMSs in the spectral type range K0 - M5. These 



are the main results: (1) the fraction of stars with ongoing mass 
accretion decreases fast with time, going from ~ 60% at 1.5-2 
Myr down to ~ 2% at 10 Myr; (2) this fraction is systematically 
lower than the fraction of stars showing near-to-mid infrared ex- 
cess; (3) mass accretion and dust dissipation in the inner disk ap- 
pear to have different characteristic timescales, 2.3 and 2.9 Myr 
respectively; (4) within 5 Myr the mass accretion rate of 95% of 
the stellar population drops below our detection limit of 10" 11 
M G yr~ l . While viscous evolution and photoevaporation might 
be unable to slow down accretion (and leave a substantial dust 
mass) on such a short timescale, planet formation, and/or migra- 
tion, in the inner disk (few AU) might be a viable mechanism to 
halt further accretion onto the central star. 
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